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Abstract: We study the process of tachyonic preheating using approximative quantum equations 
of motion derived from the 2PI effective action. The 0(N) scalar (Higgs) field is assumed to 
experience a fast quench which is represented by an instantaneous flip of the sign of the mass 
parameter. The equations of motion are solved numerically on the lattice, and the Hartree and 1/N- 
NLO approximations are compared to the classical approximation. Classical dynamics is expected 
to be valid, since the occupation numbers can rise to large values during tachyonic preheating. We 
find that the classical approximation performs excellently at short and intermediate times, even 
for couplings in the larger region currently allowed for the SM Higgs. This is reassuring, since all 
previous numerical studies of tachyonic preheating and baryogenesis during tachyonic preheating 
have used classical dynamics. We also compare different initializations for the classical simulations. 
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1. Introduction 

In models of hybrid inflation (l], 0j, one scalar field (the inflaton, a) triggers the symmetry breaking 
of a second (the Higgs field, 0) 1 at the end of or after inflation || [|. Typically the inflaton-Higgs 
coupling is introduced as an effective contribution to the quadratic Higgs term, 

\^4> 2 = \{~^~^ 2 )4> 2 - (hi) 

As the inflaton goes down, this effective mass nl s (t) changes sign from positive to negative, and as 
a result the Higgs field performs symmetry breaking by "rolling" into its broken phase minimum. 
The rolling off its unstable maximum is accompanied by an exponential growth of the amplitude of 
the momentum modes of the Higgs with fJ^ s (t) + k 2 < 0, a process known as "spinodal instability" 
or "tachyonic preheating" Interactions end the transition, and the energy becomes redistributed 
through scattering into all the modes and other fields coupled to the Higgs. Eventually the system 
will thermalize to some temperature, which is called the reheating temperature of the Universe after 
inflation. Symmetry breaking and preheating are nonperturbative out-of-equilibrium processes, 
which are difficult to treat analytically {| . 

1 Although our main emphasis is on the Standard Model Higgs, the discussion applies to a general scalar field. 
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Numerical studies of tachyonic preheating || ||, [lC , 11, 12| and baryogenesis during tachyonic 
preheating jl3|, [l4| apply the classical approximation to the out-of-equilibrium field dynamics, in 
which the evolution is determined by the classical equations of motion. The argument for making 
this approximation is that the important unstable low momentum modes acquire large occupation 
numbers; see |ll| [ll| for a detailed justification. 

Whereas the instability and growth of the scalar field modes is due to the mass squared [i^g 
becoming negative, it is not obvious what happens when the Higgs is coupled to gauge fields. In 
O, [b| it was seen that the low-momentum modes of the gauge field also acquire large occupation 
numbers, supporting the classical approximation for the coupled gauge-Higgs system. This large 
growth is essential for successful baryogenesis during (electroweak) tachyonic preheating |l7|, [ll| ||, 
|l3| , [l4| , since these low-momentum modes are supposed to dominate the dynamics of the anomalous 
baryon number violating processes in the Standard Model (SM) (for a review, see Jl|] ) . 

In an attempt to use a quantum description, while still incorporating the inhomogeneous nature 
of the typical non-perturbative field configurations of defects and anomalous baryon number vio- 
lations, fermions in classical backgrounds were studied i n ]2C|, pl| , and scalar fields in the Hartree 
approximation with inhomogeneous mean fields |22| |24j, B5j, in 1+1 dimensions. 

In this paper we apply an alternative description of out-of-equilibrium dynamics of quantum 
fields, using 2PI effective action techniques |2(| ^7|. In this approach equations of motion for 
quantum correlators are derived from a suitable quantum effective action (the 2PI effective action) , 
in which the basic variables are the one- and two-point functions. Although the approach is formally 
exact, practical calculations require us to truncate an infinite expansion of loop diagrams ($), 
defining a truncated 2PI effective action. The stationarity of the truncated effective action provides 
the <f>-derived equations of motion. These equations include infinite loop resummations in the same 
manner as truncations of the Schwinger-Dyson hierarchy, and the two are closely related p8| (for 
other related work see @ |(| ^> H> H)- 

A crucial feature of ^-derivable approximations is that the derived equations of motion conserve 
global symmetries and associated Noether currents in time. In particular the energy is conserved. 
This is a useful feature when studying out-of-equilibrium processes where most other quantities 
such as particle distribution functions evolve in a complicated way. The conservation of energy 
applies not only to the full effective action, but to any level of truncation. 

In p7{ the formalism was numerically implemented for a relativistic scalar field in 1+1 di- 
mensions for a study of the approach to equilibrium. Subsequent studies include |j4| |35| ] in 1+1 
dimensions, |}6| in 2+1 and (3^, |3^| for fermions and scalars in 3+1 dimensions. The quantum 
<I>-derivable scheme was compared to full and ^-derived classical approximations in J39|. Finally, 
the process of resonant preheating was studied in [[hi] . 

We report here on a study of tachyonic preheating in an 0(4) model of scalar fields similar to 
the Higgs sector of the SM. However, since there are no gauge fields present, we expect Goldstone 
modes to appear. In the SM these would turn into the longitudinal modes of the massive gauge 
bosons. The results presented here may also have a bearing on effective pion dynamics in heavy- ion 
collisions. 

In particular we will be interested in comparing results from the ^-derived equations to those 
obtained in the classical approximation. Initially we tried truncations of $ based on the weak 
coupling expansion, but the results were problematic because we run into instabilities. We record 
here results based on an expansion in l/N, with N the number of fields (here N = 4), following 
|35| , Ell. This expansion is systematic and the numerical algorithm is stable even for very large 
occupation numbers |55|, [hj| . 
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2. Tachyonic preheating 



The system is described by the classical action 



S 



d 4 x 



— (j3a{x)(j)a{x) + 



(Mx)Mx)) +V a 



(2.1) 



with Vq = 3A^/i 4 /2A the energy density at (f> = 0, and a — 1, ••• ,N, summation implied. In 
equilibrium at zero temperature, the system is in a broken-symmetry phase if fj, 2 > 0, with 



A 



(2.2) 



where n a is a unit vector, n 2 = 1. The particles in this phase are N — 1 massless Goldstone bosons 
corresponding to the modes transverse to n a , and the Higgs boson corresponding to the longitudinal 
mode, with mass m = V2 /i. A feeling for the strength of the interaction may be obtained from the 
amplitude for the scattering of Goldstone bosons a + b — ► c + d (a,...,d= 1, . . . , N — 1, assuming 



A (Pa +Pb) 2 

3iV (Pa+Pfc) 2 +™ 2 



2 perm. = S ab 6 cd 



(Pa+Pb) 2 



0(P 4 ) 



2 perm. 



(2.3) 



This leads to an S-wave unitarity bound A < 48irN/(N + 1) 



For N = 4, the triviality bound 



obtained with the lattice regularization is about two-thirds of this, A < 80 J42, 43 1. In the Standard 
Model, v = 246 GeV, and a Higgs mass m = 175 GeV would give A = 6. This looks like a weak 
coupling, but such a Higgs mass is near the upper end of the interval set by the constraints from 
radiative corrections |Q. 2 

The coupling of the model to the inflaton may be described by the replacement — fi 2 — > fi 2 s in 
the action (|2.l| ). The spinodal instability occurs when the coefficient of the quadratic term in the 
potential flips sign. Ref. |i"5| contains a detailed analysis of the initial subsequent development, and 
a justification is given of a classical approximation after the instability has sufficiently progressed; 
see also p5[ . We go to the limit where the flip happens infinitely fast, a quench, which we model 
by 



Mcff(* < 0) = f 



> 0) = —fJ 2 - 



(2.4) 



We assume that the initial state of the system is the ground state in a quadratic potential just 
before the quench, a free "vacuum" |0)o, with correlators 3 



(0K<^ k |0)o 



(OKttMOo 



2sjy? + k 2 



Sab, 



VM 2 + k 1 



Sab 



(O|0k^-k|O) O =Jab, 



(2.5) 

(2.6) 
(2.7) 



where k = |k|. As detailed in pq| , we can then solve for the evolution of these correlators in 
the quadratic approximation (A = 0) as a function of time. The result is reproduced here for 



2 On the other hand, a potential application to heavy-ion collisions using effective pion dynamics, where the role 
of the Higgs is played by the sigma-resonance, m — > m CT « 600 MeV, v — > = 93 MeV, would require very strong 
coupling, A a; 120. 

3 We use the Fourier decomposition <f>{t, x) = L~ s ' 2 5TJ k 0k(*) e!kx m a periodic box of linear size L. 
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convenience, 



cr(t)s ab =(o\K(tW- k (t)\o)o = 



_ n k {t) + 1/2 



OJ k (t) 



, + 2 



Cr(t)S ab =<0|7rg(t)7ri k (t)|0)o = [n k (t) + 1/2] Wfc (t) 5 ab 



,+2 



1+ co s W) 



<5afc, 



Cf*(i)<U =<O|0£(t)7ri k (t)|O) o = [n*(t) + i/2] 5 a6 



, +2 



1 sin(2w fe i) + 



Sab, 



(2.8) 



(2.9) 



(2.10) 



with ujfc = \J ±/i 2 + k 2 . The derived occupation numbers and frequencies are 



n fc (t) + l/2=^Cr(t)Cf(i), afc(i) = ^C^(t)/Cf(t), ri k (t)=ReCt(t). (2.11) 



It follows that modes with /c 2 < /x 2 grow exponentially with 



njt(f) oc exp (2a/ /i 2 — k 2 t). 



(2.12) 



We will use equations (2.8) - ( 2.10p for comparison at early times. 

Also in more general situations (2^, [l6| ^ |l2| the above u) k give an instantaneous charac- 
terization of frequency assuming a translation and rotation-invariant density matrix, and they 
provide for a corresponding definition of annihilation operators as (concentrating on one real field) 
<2k = ^= (wfc^k + iTTk), and creation operators <z£, such that (a^Ok) = n&. Hence the n k have a 
natural interpretation as occupation numbers. For free fields in equilibrium the n k and Lo k coincide 
with the standard particle numbers and energies. The h k are given by Im(aka-k) = n k . For 
large n k the non-zero commutator corresponding to the imaginary part of ( 2.10| ) becomes unim- 
portant in expectation values of suitable observables (e.g. particle numbers), which suggests that 
the classical approximation should be good. Ref. jy| stresses that F k = n k 1 is the important 
criterion for classical behavior. For our initial state n k and h k are not independent, but satisfy 
(n k + 1/2 + fi k ){n k + 1/2 — n k ) — 1/4, so with h k — \/n 2 +n k , n k > \/2 whenever n k > 1, and 
for large n k the two become very close. For precise formulations of the classical approximation in 
terms of probability distributions of initial conditions, see |lj| jig ]. 

The initial state is O(N) symmetric, which implies that the mean field stays equal to zero also 
for the general case A / 0, 



(f> a (x) = 0, a =1, ■••,JV, 
and that the correlators can be written as 

{T(j) a {x)(j) b {y)) = 6 ab G(x,y). 



(2.13) 



(2.14) 



As the n k grow, eventually the back-reaction from the quartic term will become important and the 
evolution will deviate from the quadratic approximation. As we will see below (section (0)), this 
happens around the time when (t = x°) 



M 2 {t) = - f i 2 + X 



n + : 

6N 



[G(x,x)-G(0,0)}~0. 



(2.15) 
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Since 

G(x,x)-G(0,0)c f^^ y (2.16) 

we expect the occupation numbers to be non-perturbatively large 

n k ~ 1/A, (2.17) 

which suggests that a simple loop expansion in powers of the coupling is unreliable for tachyonic 
preheating. 

3. Broken phase and finite volume 

After the initial tachyonic preheating the system will go to equilibrium and thermalize. To estimate 
the final temperature T, consider a simple approximation in which the initial energy density, Vo , is 
transfered to a free gas of N — 1 massless Goldstone bosons and one massive radial mode of mass 
77i = y2fj,. This approximation is expected to be reasonable for weak coupling, provided T is far 
from the critical temperature T c at which the transition between the broken and symmetric phase 
occurs. Using Vo — 3N{i 4 /2A, the temperature T is then given by 



3iYy f d 3 k (N- l)k v/2/i 2 + k 2 



2A J (2tt) 3 \ e k / T - 1 eV /2^+kyT _ 1 



(3.1) 



For N — 4, A = 6 this gives T — 0.95 fi (the large N limit would give 0.93 /u). On the other hand, 
the transition temperature T c can be estimated in a simple one-loop approximation |47| , 

367V ^ 1/2 



T ^{WT^) m ' (3 ' 2) 

which gives T c = 2.83 fi, for N = 4 and A = 6. This also agrees well with Monte Carlo results |47| . 
Since the estimated final T <C T c , we expect the system to end up in the broken-symmetry phase. 

Our simulations are performed in a periodic box of size I? . It is well-known that quantum 
tunneling effects prevent spontaneous symmetry breaking in finite volume. Hence, (<p a ) = 0, but 
this does not mean that physical effects of symmetry breaking cannot manifest themselves in finite 
volume. In fact, the Goldstone bosons may cause relatively strong finite size effects and slow 
thermalization. 

The standard way to define (<p a ) is to add a symmetry breaking term to the action, 



AS = J d xen a <fi a , (3.3) 

which pulls (4> a ) into the direction n a , and then evaluate 

vn a = lim lim (<f> a ). (3-4) 

e— >0 L— >oo 

Aided by finite-size scaling analysis this method has been used fruitfully in Monte Carlo studies 



Jlq ]. The order of the limits in (3.4) is important. However, reversing the order (working at e = 0), 
a lot can be learned from analysis at finite size [ |49"| . 

A proper finite-size scaling analysis is outside the scope of our present explorative study. In- 
stead, we shall work at reasonably large volumes, L\i > 10. As a guideline for interpretation of 
the numerical results at relatively large times (close to equilibrium), we take the tunneling into 
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account by simply averaging the infinite- volume broken-symmetry correlators over all internal 0(4) 
directions: 



(0a) = vn a -> 0, 
(T(f>a{x)(t>b{y)) = v 2 n a n b + (S ab - n a n b )G G (x,y) + n a n b G a (x,y) 
-> 5 ab G(x,y), 

"' 2 N - 1 ~ , s 1 



G(x,y) = - 



y ■G G (x,y) + -G K (x,y), 



(3.5) 

(3.6) 
(3.7) 



where G and H denote the Goldstone and Higgs contributions. Hence, after having settled sufficient 
into equilibrium, the O(N) symmetric Green function G is expected to have both Goldstone and 
Higgs contributions, and a zero-momentum mode expressing the condensate. 

4. The 2PI effective action and equations of motion 

The 2PI effective action can be written in the form [E6| : 



m G] = S[4>] - l -Tr InG + *-Tr G^G + $[0, G], 



(4-1) 



with 



xJr TT1 \ = ( d * dab + V 2 Sab - 7^7 ^ c Sab ~ ^7 ^ b ) 6 ( x ~ y)- ' 42 ! 
5<f> a {x) <p b {y) V 3iV ' 



We use the Schwinger-Keldysh formalism, in which the fields live on a contour C running from 
t = to some time t max and then back to t = 0. Time integrations and time ordering are along 
this contour. At the two times t± — there are functional integrations implementing the density 
matrix that specifies the initial state. The functional G] is a sum of 2PI skeleton diagrams 
with bare vertices and full propagators corresponding to the initial state , see figure [l]. 



i <*> = 



Figure 1: Expansion of $ in terms of diagrams. The lines denote the full G, the vertices denote the bare 
vertex functions iS n S[(j>]/5(j>i ■ ■ ■ 6<j> n , n = 3,4. Symmetry factors are not indicated. 

The full propagator is a stationary point of the effective action. Taking the functional derivative 
gives the equation determining G ab , 






5G ab (x,y) 
with £ the self-energy, 



= °^ G ab( X >y) = G Q, 1 a b ( X ^y)+ i ^ab(x,y), 



H ab {x 7 y) = -2 



8G ab (x,y)' 



Multiplying ( |4.3| ) by G (y,z) turns it into an integro-diffcrential equation, 



/ dt' 


[d 3 y 


IC o 





ibe / 



(4.3) 



(4.4) 



(4.5) 
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Together with the equation for (j>, which follows from 6T[<j>, G]/5cf) a — 0, this determines G(y,z) 



in a self-consistent way, since the self-energy is itself a functional of G, eq. (|1J). Since will be 
identically zero because of the O(N) symmetry of the initial state, we now specialize to this case. 
The delta function in ( |4.5| ) for G can be taken care of by introducing G > and G < according to 

G ab (x, y) = 6 ab [G> (x, y) 6 c {t-t>) + G< (x, y) 6 C (t> - t)] , (4.6) 

where the step functions 9c (t — t') refer to the time label along the contour C. The 'statistical 
function' F and the spectral function p correspond to the real and imaginary parts of G > or G < , 

G>(x,y) = F(x,y)-ip(x,y)/2, (4.7) 
G<(x,y)=F(x,y)+ip(x,y)/2. (4.8) 

In operator language they correspond to 

F(x,y)S ab = ^({r(x)^ b (y)}), p(x,y)S ab = i<[^-(ar), 6 (y)]>. (4.9) 
The functions F(x,y) and p{x,y) are real and have the symmetry properties 

F(x,y) = F(y,x), p(x, y) = -p(y, x). (4.10) 
The canonical commutation relations imply 

dop{x,x')\ x o= x ,a = -dwp(x,x')\ x o= x ,a = £(x,x'). (4-11) 



Equation (4.5) now leads to coupled integrodifferential equations for F and p, supplemented by 



initial conditions depending on the initial density matrix. These equations preserve the properties 



eq. ( 4.10 , 4.11 ), and conserve energy for any choice of truncation of $[</>, G]. 

One would like to know if truncations of <& lead to solutions for G that describe massless 
Goldstone bosons when the system is in the broken-symmetry phase . Generically this is not the 
case 

However, one can introduce the 1PI effective action 

f[0] = r[0,G[0]], (4.12) 



with G[<f>] the solution of (4.3) for given mean field <fi, and define a Green function G (called 'external' 



501) in the usual way from T: 



^'■d- ^MWrt - (4 - 13) 



Then the usual derivation of Goldstone's theorem in terms of V applies to G _1 , it has N — 1 zero 
eigenvalues at zero momentum p| , |50| |. Without approximation the Green functions G and G 
coincide, but this is no longer true for a truncated <&. The difference between G and G depends on 



the functional derivative 6G[(f>}/5<fi |50 . This is a three-point function (given by the solution of a 
non-trivial integral equation) that vanishes in our case of zero mean field. Therefore, we do expect 
G to describe massless Goldstone bosons. 

At first we tried a truncation of the effective action based on the loop expansion, by keeping 
the diagrams shown in figure [j], for N = 1. However, this lead to severe instabilities before the 
transition was completed, for weak and also for stronger couplings, which we were unable to cure. 
We interpret this as a breakdown of an expansion based on the smallness of A, because of the 



non-perturbatively large occupation numbers expected during the instability (eq. ( 2.17 )). We then 
turned to an expansion of the 2PI effective action, not in powers of the coupling A, but in terms of 
1/N. Following Q we expand $ to Next-to- Leading Order (NLO) in 1/N. The explicit form of 
the resulting equations can be found in pj| (see E| for the case of non-zero mean field). 

4 Such Goldstone bosons are still present in our case of zero mean field, cf. section ^[ 
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5. Numerical implementation 



We discretize the O(N) model on a space-time lattice with the action 



Siat = a s «t 2^ I 2 2 2"^a(^ x ) ~ 247v (^a^.x)) (5.1) 



Here the lattice spacings are a s and at in the spatial and time directions, dt4>(t,x) — (<fi(t + 
at, x) — x))/a t , and similarly for spatial derivatives. We use periodic boundary conditions, with 
volume is L 3 = (n s a s ) 3 . In the following we shall use lattice units, a s = 1 and write di for the 
dimensionless time-step (dt = at/a s ). The coupling Ao and mass parameter Pq are bare parameters 
to be determined below. At tree level, Ao = A and p^ = —p 2 in the broken phase. The lattice 
version of the squared momentum is given by, 

3 o 

= ^(2-2 cos fc,), h = 1, ni = -^- + !,-■■ 5.2 

z — ' n s I 2 

i— 1 

(n s is chosen even). Plotting data as a function of the lattice momenta rather than continuum k 2 
corrects for most of the lattice artifacts. 

As explained in , truncating the 1 /N expansion at NLO and specializing to a homogeneous 
system (G(x, y) — G(t, t', x — y)) with zero mean field, leads to the equations of motion: 

d t d' t F(t, t 1 , x) = dt&iF{t, f, x) - M 2 (t)F(t, t', x) 
t 

t"=0 z 
{' 

-^2 dt^2E F {t,t",z)p{t",t',x-z), (5.3) 

t"=0 z 

d t d' tP {t, i',x) = did' iP {t, i',x) - M 2 (t)p(t, t', x) 
t 



with d t 4>(t,x) — (<f)(t,x.) — <j)(t — dt, x))/dt, and similar for and with 5 



M 2 (t)=p 2 + \ ^^F(t,t,0), (5.5) 
£ F (t, t", z) = ^ (F / F - pI p /4) (*, t", z), (5.6) 



X p (t,t",z)=-^L(FI p +pI F )(t,t",z). (5.7) 



5 When practical, we shall use a notation in which, for instance, p(t, t' , x)F(t, t' , x) = (pF)(t, t' , x). 
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The functions I F and I p represent the summation of bubble diagrams included at NLO |35| 
/ F (M',x)=^(F 2 -p 2 /4) (*,*', x) 



t"=0 

t' 



I„(t,l/,x)=-±(Fp) (t.t'.x) 



A„ 



£ dt^2l p (t,t",z) (Fp) (*",*', x-z). 



(5.8) 



(5.9) 



Note that £f and E p contain an overall factor A 2 ,, they would start at order A 2 , upon solving the 
above equations by perturbation in Ao. The contribution of order Ao (coming from the first diagram 
in figure Q) is purely local (non-zero only for t = t") and contained in M 2 defined in (5.5). During 
the numerical evolution we have to solve for I F and I p parallel to solving for F and p. Although eqs. 
(5.8, 5.9) look implicit, by doing things in the right order and taking advantage of the symmetries 

I F (t,t',x)=I F (t',t,x), I p (t,t',x) = -I p (t',t,x), I p (t,t,x)=0, (5.10) 

they can be solved explicitly. The energy functional is given by 

p* (( , 0)) . + &M|M^ + & F «, t, o) + FHt.,., o) 



A 



12 N 



E dt J2 {2pFI F + (F 2 -p 2 /A)I p ) (t,t",z)-sub., 



(5.11) 



where 'sub.' is a subtraction. Ideally, this subtraction is such that the energy density is zero in the 
vacuum (the zero-temperature ground state), but its evaluation requires a non-trivial minimization. 
In practise we set the energy density to Vq, for the initial state (pTs| 2.7). 



5.1 Hartree approximation 

In addition to the NLO-l/N approximation, we study the Hartree approximation, which results from 
taking into account only the first diagram in figure 0. This means leaving out all time integrals, 
and the I F , P , keeping only the first line of eqs. ( |5.3| ) and (|5.4D, with M 2 defined by eq. (|5.5|). Notice 
that this is not the Leading-Order (LO) approximation in 1/N, since the coefficient of the local self- 
energy is Aq-^2 as opposed to ^ for LO-l/N. Neither of the two includes non-trivial scattering 
in the dynamics and they give qualitatively the same results. We use Hartree dynamics only for 
comparison. 



5.2 Initial conditions 

We restrict ourselves to initial conditions for which the density matrix is gaussian. This assumption 
is equivalent to the initial state being completely described in terms of (one- and) two-point functions 
]2j| . Consequently, an initial condition is a choice of initial F's, p's and (which in this case is 
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zero). We need to specify 



i({0 Q (O, x), 4>\Q, 0)}) = Sab F(t, t', x)| t=t , =0 , (5.12) 
i({^(0,x),^(0,O)}) = S ab d t F(t,t\x)\ t = t ,= , (5.13) 
i({7r a (0,x),^(0,O)}) = 6 ab d t d t ,F(t,t , ,x)\ t=t , =() , (5.14) 



in addition to 



i([0 a (O,x),^(O,x)]) = W(0,0,x) =0, (5.15) 
z([7r a (0, x), /(0, 0)]) = S ab 8 t p(t, t', x)| t=t / =0 = <U <5 x ,o, (5.16) 

(5.17) 

where we used G(x,y) = G(i,£',x — y) for a homogeneous system. This also allows us to go to 
momentum space and define 6 

1 ^init i i In 

5<{«(0),^ k (0)}> = <UF k (0,0) = k XJ- S ab , (5.18) 

^(K(O),0 b - k (O)}) =6 ab d t F k (t,t%= t , =0 = h{ nit S ab , (5.19) 

^({<(0),^ k (0)}> = <UdAF k (M')|t=t'=o = K nit + l/2)4 nit 5 afc , (5.20) 

given a choice of "occupation number distributions" nP lt , n™ lt and dispersion relation w™ 1 *. We 
specialize to the case n™ 11 = 0, noting that in equilibrium this correlator is zero. As mentioned 
in section 0, we assume that the system is initially in the state corresponding to the ground state 
before the quench, which we approximate by the free-field vacuum in the symmetric-phase, 

4 nit = 0, o^^ + fc^. (5.21) 

Notice that we use the renormalizcd quantity /j 2 in the initial condition, in accordance with section 
|^. Also, since this initial state is a free-field vacuum, it is not the true vacuum of our interacting 
system in the symmetric phase. 

5.3 Choice of bare parameters 

The renormalizability of ^-derivable approximations was studied in |5l], |52|, [54|, |55| . We assume 
the NLO-l/N approximation to be renormalizable, in the sense that we can choose the coefficients 



in the discretized equations (5.3-5.9) in such a way that the lattice spacing becomes exceedingly 
small 7 - loosely called 'the continuum limit'. 8 We will not study lattice-spacing dependence but 
do want to choose the parameters /ig and Ao such that the relevant length scale in our simulation 
is larger than the discretization scale, i.e. /i -1 > a s . 

For this purpose we implement renormalization at one loop order. At weak coupling, mass 
renormalization is important, but not coupling renormalization. Two-loop effects are also small, 
and the bare mass ^Lq in the equation of motion is therefore estimated in one-loop renormalizcd 



6 The Fourier transforms are given by F k (t,t') = L -3 / 2 F(t, f', x) e _lkx , etc. 

7 This may imply a more elaborate set of parameters than just /j 2 and Ao- At worst, an independent parameter 
may be needed for every contribution to the truncated <J> (e.g. every diagram), since there is no symmetry relating 
various contributions. In practise this is not a problem for sufficiently accurate truncations [ ps| . 

8 Such terminology ignores triviality. 
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perturbation theory as 
ci(a s fi) = -i 



N ■ 



6N 

T/as d 3 k 



2 f \ 1 



= C + C 2 a 2 /i 2 + 



dko 
~2~7 



-2 1.2 



lat 



-* (2^) 3 2 ^M 2 + A£ 



lat 



16tt 2 



s Vln(a 2 Ai 2 ) + 0(a^ 4 ), 



(5.22) 
(5.23) 
(5.24) 



in the continuous-time and infinite- volume limit. In ( |5.22 ) the plus (minus) sign refers to the 
symmetric (broken) phase. The one-loop self-energy refers to the symmetric phase, the broken 
phase should give the same divergent contribution as a s — > and a small difference in the finite 
contribution. An idea of the size of the coefficients may be obtained from a calculation on a 
symmetric (at = a s ) lattice in imaginary time pq|, 



Co = 0.155. 



C 2 = 0.0303 . . 



Similarly, the bare coupling can be estimated by the one-loop expression 

,N + 8 / I , , , os 1 



A n = A - A 



67V 



16P 



In(afy) 



16P 



-a 



(5.25) 



(5.26) 



(where A is defined as the value of the four-point vertex function at vanishing external momenta). 
For the maximum coupling that we use (A = 6) and lattice spacing 0.5 < a s [i < 1, the difference 
between Ao and A is less than 10%. In practise we simply choose Ao as if it were the renormalized 
coupling. In our simulations we use the bare coupling in ( 5.22 ), i.e. on a n^-site lattice, in lattice 
units, 



= - - A 



N + 2 

6N 
1 



ci(^), 



hat 



0.21. 



(5.27) 
(5.28) 



Given the input parameters fj^ and Ao, the output physics is of course not known precisely, it is 
determined by the <I>-derived equations. For example, the mass of the (unstable) Higgs particle will 
differ from the tree graph value m = \/2 /i. Similarly, the renormalized coupling will not satisfy the 
relation ( 5. 26] ) precisely (A may also be defined in terms of the vev by v 2 = 37Vm 2 /A where the vev 



may be identified from the low-energy scattering amplitude of the Goldstone bosons, cf. (2J3)). But 
the differences are expected to be small, and they can in principle be determined numerically. 

In our out-of-equilibrium study there is another dimension-full scale: the mass fi in the initial 



conditions (5.21). It is the renormalized mass in the symmetric phase, and it is this mass that we 
use in (5.27) for choosing /ig. We have to keep in mind that the gaussian approximation in the 
initial conditions and the suddenness implied by the quench may lead to divergencies in the pressure 
near t = 0, in the continuum limit, see e.g. j57|. These may be dealt with by adjusting the initial 
density matrix, but it would be better to avoid the instantaneous quench and use a more realistic 
model of the transition. Such divergencies do not appear to be important in our simulations at 
finite lattice spacing. 

5.4 Observables 



We follow the evolution of F(t, t, 0) and the energy, eq. ( p.ll| ), and calculate occupation numbers 



rife and frequencies Uk through the inverted version of eq. (5.1 



n k + 1/2 =c fe VftcVF k (M')|t=t'^k(M), 
u k =Vd t d t <F k {t,t% =tl /F ]i (t,t). 



(5.29) 
(5.30) 
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The quantities and ujk coincide with the true particle occupation numbers and frequencies when 
applied to a free field in equilibrium, and have proven to be very useful out-of-equilibrium in 
interacting theories as well [l^] . The correction ct takes care of errors associated with 



the time discretization on the lattice. It is given by 24 



Ck 



(5.31) 



In the present case we argued in section |^ that the two-point function G receives contributions 
from both the Higgs and the Goldstone modes, so in this respect the and u>k are compound 
observables. 

After some initial transients, the system will approach thermal equilibrium, although this may 
take a long time because of critical slowing down: the scattering of the Goldstone bosons becomes 
suppressed for momenta smaller than the Higgs mass, cf. (^^). The approximate expression (3.7) 



for G in terms of G H and G G provides a guide to what form we may expect for the rife and u>k 

N-l n£(t) + l/2 , 1 <(t) + l/2 



defined in (5.29) 



v 2 (t)L* 
Fk{t,t) — — d k: o 



d t d t ,F k (t,t')\ t=t 



N 
N - 1 



N 



Y ("?(*) + 1/2) ^(t) + - (n£(f) + 1/2) w H(t), 



(5.32) 
(5.33) 



where, assuming that the individual Higgs and Goldstone modes have thermalized to a Bose-Einstein 
(BE) distribution, 



1 



/T_i 



e^/ T - 1 ' 



(5.34) 



and v and m are to be interpreted as the condensate and effective Higgs mass at temperature T. We 



have allowed for a time dependent v(t), although in ( 5.33 ) we have neglected its time derivatives. 
Close to equilibrium we expect the change in v to be very slow (v/v <C fJ., see figure ^ below). 
At low momenta (fc < to) the Goldstone bosons will clearly dominate. Note that in finite volume 
the condensate appears only in the zero mode of F, and there is no zero mode in the Goldstone 



contribution G . it is represented by the global rotations making up the 0(4) average in (3.6) 





Figure 2: Left: Compound frequency spectrum from ( 5.2£ )— (5.34), uj k versus k, for N = 4, T = 0.95^, 
m — v2[i, in units fx = 1. The dashed lines u> 2 — 1.395 k 2 and u> 2 — 0.5/i 2 + k 2 are the asymptotes (5.35) 



and (5.36). Right: Compound particle number, ln(l + 1/n.fc) versus Wfc. The dashed line represents u>/T. 



In figure [I (left) we have plotted the resulting u>k at the expected temperature T = 0.95 /i found 
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in section || for N = 4, A = 6. The straight lines correspond to the asymptotes 



m(n? + 1/2) 

r(iv-i) 
2 + ^ + o(fc- 2 ), 



fc 2 + 0(fc 4 ), 



fc 2 > m 2 ,T 2 . 



(5.35) 
(5.36) 



This latter formula may provide an estimate of the effective Higgs mass, although the emphasis on 
the large momentum region is questionable. A better way to determine m should be to find the 
position of the peaks in the Fourier transform of the spectral function, as done in 1+1 dimensions 
in |3. Figure || (right) shows the compound particle number from (|5.29 )-( 5.34 ) for the same 
parameters as for Wk in the left-figure. Plotted is ln(l + 1/fifc) versus LOk, which would be a straight 
line with slope 1/T if rik were a Bose-Einstein distribution with zero chemical potential. The 
deviation from BE is hardly visible, which suggests that such a plot is also useful for judging the 
approach to equilibrium, even for the compound particle number. For a BE with a non-zero chemical 
potential \i ci we would see an intercept with the w^-axis equal to fi c . For the compound particle 
number this also applies with the understanding that the intercept gives the chemical potential of 
the Goldstone particles, /i^. 

For sufficiently large volume the condensate follows from the zero mode of F, 



v 

N 



F (t,t) 
L 3 ' 



(5.37) 



which has been found to give good results compared to a more rigorous procedure based on (3.4) 
and finite-size scaling [jt8[. The relative correction 



+ 1/2 \(n$ + 1/2) 



>J 2 L 3 6\/2A^L) E 



(5.38) 



(which may also be estimated from dtdt> Fo(t,t') t= t' / Fo(t,t)m 2 ) is tiny: for /iL = 10 and An[/ of 
order 100 it is only of order 10~ 2 /N. 

In the classical approximation we found |lq, 14, 12| the observable N~ l {4> a {x)4> a ix)) useful 



for monitoring the tachyonic transition, and we have used it here as well. Its quantum analog 
F(x, x) = F(t, t, 0) is divergent in the continuum limit and needs a subtraction. We have taken this 
to be its initial value F(0,0) = F(0,0,O) = ci(^) (~ 0.21 in lattice units cf. fTl^), (|5~2ll ), (^28|)). 
The difference F(x,x) — ^(0,0) then vanishes at time zero and at late times ( 5.32 ) indicates that 



F(*,*)-F(0,0) = ^ + i,W 

3 k \ 



N- 1 nf 



1/2 



v 

iV 



N 

N ■ 



1 n\ 



(1 



1 ng + 1/2 
^ 0) + N 




(2tt) 3 
d 3 k 



N 



N LO? 



N - 



7r /a. 



(2tt) 3 



N 



1 1 

hat 



N V™ 2 + kL 



^lat , 



(5.39) 



where the second line refers to a s ^ 1 units. The last integral is logarithmically divergent as a s — > 
(at two-loop accuracy there may be also quadratic divergencies). A different subtraction, e.g. such 
that the result is finite in the broken-symmetry vacuum state, may render it divergent as t — > 0. 
This is an example of the artifacts introduced by the suddenness of the quench. Nevertheless, for 
the lattice spacings that we use these remaining divergencies are negligible. For example, in lattice 
units the last integral in ( 5.39 ) is only 0.00922 for m 2 = 2/i 2 



1. 
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We also compute the 'memory kernels' in the equations of motion ( p.3[ )-(5.4), the self-energies 
~^Fk(t,t') and E p k(i,i')j which can be compared with perturbative estimates. For long runs, finite 
computer memory requires us to "cut off" the self-energy (in terms of the functions If,p), by keeping 
the memory kernels only some finite range backwards in time (i.e. £(t, t",z) — > 0, \t — t"\ > t cu t)- 
The size of the self-energies will help us determine whether we have cut off late enough for the 
discarded memory integrals to be negligible for the dynamics. A too early cut-off also shows up as 
a non-conservation of energy. In all the runs presented here, the energy is conserved to within two 
percent. At early times its fluctuations are due to time discretization errors. At very late times, 
the effect of cutting off the memory kernel can add up to a very slow drift in the energy, but still 
within the two percent band. Keeping the whole kernel (say, for lattices of smaller spatial extend) 
there is no such drift, and we have checked that later cut-offs make the drift smaller. Our choices 
of memory cut-offs reflect a compromise between exact energy conservation, lattice sizes, computer 
memory and especially CPU time, and the wish to study reasonably long time evolution. 

6. The classical approximation 

In the classical approximation, a set of classical initial realizations is evolved using the classical 
equations of motion. These realizations are drawn from an ensemble of initial conditions that 
reproduce the initial quantum correlators in a classical approximation. In the case of tachyonic 
preheating we do not have to wait for the correlators to become classical if the couplings are weak 



we can go back all the way to time zero and reproduce eqs. (2.S), (2.9), ( p.lCj ) at t = 0, 



i , / ..i T p 

Ct*(p) = — ^4 , CT(0) = V o , k<fi. (6.1) 

See fni for a more detailed justification of this procedure. Note that the stable modes k > [i 



are initialized to zero. In section 7.3 we comment on the case where all the modes are initialized 
according to the initial quantum correlators, including the modes k > fj,. Observables are then 
averaged over the sample, giving an approximation of the quantum observable, 

(O) 

quant. ^ (0)class.cn SC mblc (6.2) 

To compare with the <I> derived approximation, we perform classical simulations starting from 



the same lattice Lagrangian as in the quantum case, eq. (5.1). The classical equation of motion 
reads 

a&0„(t,x) ^V^iftxj + AM) - A (0 b 6 a ) (i,x). (6.3) 
^-^ 6 N 

The classical simulations use the same lattice sizes and spacings as in the quantum case. As 
observables we use 

(#t, x)^(i, x)) c i ass = / Mt> *)Mt, x) \ (6.4) 

\ S x / class 

which is the analog of the quantum correlator -F(£,x;t,x) = F(t,t,0), and the energy 



(6.5) 



Comparing (<j)(t, x)(f>(t, x)) c i ass to F(t,x;t,x) we should keep in mind that the latter is divergent in 
the continuum limit, whereas classically 

(0(0, x)0(O, x)) class = ^ Y\ ~ ; 9 X - °- 0034 - ( 6 - 6 ) 
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We therefore compare (cf. the discussion in the previous section) 



{(/>(t, x)<l>(t, x)) class F(t, x; t, x) - F(0, x; 0, x). 



(6.7) 



ignoring the small difference (0(0, x)0(O, x)) c i ass ~ 0.0034 at time zero. 

We also calculate occupation numbers and frequencies from the < 
manner similar to eqs. ( 5.29 5.30| ) (without the '1/2'), 



and 7T7T correlator in a 



n k = N~ 



LUk = 



c)/(^ 6 



(6.8) 



We average over k with the same length |k|. This gives better statistics, and since the system is 
homogeneous and isotropic as for the quantum case, the correlators depend on the length of k, 
apart from lattice artifacts. 

In the classical approximation we expect equipartition to occur after very long times, for which 
the occupation numbers for the Goldstonc and Higgs modes are given by 



T, 



class 



r. 



ha 



II f 



class 



in - 



(6.9) 



The Rayleigh- Jeans temperature is determined by classical equipartition, E = Nn 3 T c \ ass , or with 
E ~ VqL 3 , T c i ass = (3/2)(^/A)(a s ^,) 3 . In our longer-time simulation we used a s fi = 0.7 and Ao = 6, 
which gives the low value T c i ass ~ 0.086 fi. The corresponding guiding plots for the compound 
frequencies and particle numbers are shown in figure ^. The deviations from the BE case in the 
compound- frequency plot appear limited but significant. The greatest deviations appear in the 
compound-particle number case: the curvature caused by plotting ln(l + l/n/c) (instead of the more 
appropriate 1/rifc), done for diagnostic reasons, which cannot only be ascribed to a dominating 
Goldstone contribution. 

In (l+l/n) 




Figure 3: The Rayleigh- Jeans case, N = 4, T — 0.086/^, m = \/2m, units fi = 1. Left: Compound 

f 22 22222 

frequency spectru m, m k vers us k . The dashed lines ui = 3.75 k and ui = 0.5/i + k are the asymptotes 
corresponding to ( 5.35 ) and ( 5.3(: ). Right: Compound particle number, ln(l + l/nfe) versus Wfc. The dashed 
line represents the pure Goldstone case, ln(l + w/T). 



7. Numerical results 

7.1 Early times: Tachyonic instability 

We first describe simulations during the initial tachyonic transition on a lattice with n 3 = 32 3 sites, 
lattice spacing a s [i — 0.7, volume (L^j,) 3 = (22. 4) 3 , time step dt = a t /a s = 0.1, with weak coupling 
Ao = 1. The number of fields N = 4. 
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We compare four different approximations: f ) The analytical result in the quadratic approxi- 
mation (eqs. ©, Q, ( pTToj) ) ; 2) Hartree (cf. sect. 3) f/N to NLO (eqs. Q - (pUll)); 



and 4) classical (eq. ( |6.3| )). In the classical case, we average over an ensemble of 400 configurations. 
For the 1/N-NLO case relatively short times allow us to keep the whole memory kernels in the 
self-energy. 

In figure |we show the evolution of the subtracted F(t, t, 0) for the NLO, Hartree and classical 
(eq. (S.7)) case. We have also included LO-l/N for comparison. The NLO and classical agree well 
with each other (and settle later near the zero-temperature vev, v 2 /N = 6^ 2 )- The Hartree result 
is remarkably lower than the others, while the classical approximation seems to work very well. The 
Hartree being lower appears to be the result of the choice of coefficient of the local term, the choice 
of N. In the limit of large N we recover the LO result, which for N — 4 overshoots compared to 
the classical and NLO. Below, we will discard the LO approximation, since it is qualitatively the 
same as Hartree. 

Figure |^ shows the evolution of the compound occupation numbers for the Hartree approxima- 
tion (green/grey) compared to the quadratic approximation (black). At early times the agreement 
is very good but eventually back-reactions enter and the spinodal growth ends. As is well known, 
the homogeneous Hartree approximation does not include non-trivial scattering, and essentially no 
energy is re-distributed to modes of momentum k higher than /x. 

Figure H compares Hartree (green/grey) to NLO (black) dynamics. Initially they agree very 
well. As the back-reaction kicks in, scattering processes included in the higher loop diagrams of 
NLO allow re-scattering of energy into the higher momentum modes. The spinodal growth ends 
roughly around t/i = 4 — 5, when F/[i 2 ~ 4 (cf. figure^), at which time M 2 (t) defined in (5^) goes 
through zero (as anticipated in ( |2.15| )). 

In figure |?j we compare the NLO (black) result to the classical (green/grey) case. In the classical 
case only the unstable modes are initialized, which is why the stable modes are outside the picture 
initially (having nk = 0). In the quadratic approximation the evolution for classical modes and 
quantum mode functions is identical, so we expect the unstable modes to also be described by 
eq. ( 2.11 ). This is indeed what we see, and as back-reactions kick in, the stable modes become 
populated and "hook up" with the unstable ones to reproduce the low-momentum part of the 
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Figure 5: Compound occupation numbers rik vs k/fi for the Hartree case (green/grey) during tachyonic 
preheating compared to the quadratic approximation (black). The individual plots represent different times. 
Notice the logarithmic y-axis. The Hartree approximation follows nicely the quadratic behaviour through 
the instability, but back-reactions stop the growth of the modes. Very little power is scattered into modes 
with k > fi. 



distribution of the NLO case very well. Turning things around, the full dynamics of the classical 
evolution is well described by the truncated NLO quantum dynamics. 

Figure H shows the compound-dispersion relation during the transition at NLO. The low mo- 



mentum modes are again well approximated by the gaussian equations (2.8, 2.9, 2.10), and although 
the quantitative agreement is lost for large momentum modes and larger times, the qualitative pic- 
ture is well understood. For the Hartree case the picture is roughly the same at these early times. 
The averaging in the classical case was seen to smear the oscillations. 

7.2 Intermediate times: approach to equilibrium 

As interactions become important the spinodal transition ends. We may expect scattering to result 
in equilibration of the system. This does not happen in the Hartree approximation, but the inclusion 
of the NLO diagrams has been seen in previous studies to lead to thermalization, in 1+1, 2+1 with 
scalar fields and in 3+1 dimensions including fermions |^8, 3G ; 37], [58). In the present case, however, 



the suppression of the scattering of the massless Goldstone bosons at low momenta (cf. (j^) may 
lead to slowing down of the thermalization of these modes. In the classical approximation, we 
expect classical equipartition to take place after a very long time. For intermediate times, it was 
seen in Q that the low momentum modes of a system including SU(2) gauge fields could be 
well approximated by a Bose-Einstein distribution, allowing for a determination of an effective 
temperature and chemical potential, even with classical dynamics. 

To study the longer time behavior we performed simulations on a smaller lattice, = 16 3 , 
again with // = 0.7, dt = 0.1, N = 4, but with a larger coupling Aq = 6. For the classical simulations 
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Figure 6: As figure |E] comparing NLO (black) to Hartree (green/grey). At early times they agree well, but 
for NLO there is scattering into higher momentum modes k > /i. 



we averaged over 1000 realizations of the initial conditions. For the quantum NLO case we cut off 
the memory kernel at 38.5 fi . This gave no problems with energy conservation and we believe 
it does not affect the evolution much. With the enhanced back-reaction of the larger coupling we 
expect smaller n^ after the transition, since the spinodal transition does not last as long. Instead of 
a zero mode reaching previously no ~ 1000, we have here only no ~ 200. This still amply satisfies 
the criterion nu ~S> 1 for classicality. 

In figure ^ we see the evolution of the compound dispersion relation. In the NLO case, the 
oscillating behavior generated during the transition (figure ||) is washed out by scattering to a 
straight line. A close-up at the latest time the NLO compound-frequency plot figure [ll] appears to 
show the behavior anticipated in figure The asymptote to the data at larger momenta indicates 
an effective Higgs mass m 2 ~ 1.80 /x 2 (cf. ( 5. 36 )), somewhat lower than the zero-temperature value. 
At finite temperature in the broken phase we do expect a smaller Higgs mass, since the condensate 
is also smaller than at zero temperature (see figure [l^, left). Inserting the finite temperature v 2 in 
the expression for the mass m 2 = -^v 2 , one finds m 2 w 1.76 /i 2 . 

The slope of the line through the origin and the first data point is at this time 1.22, quite a bit 
off the anticipated 1.395 (see figure |J). With this slope and effective mass, no choice of T satisfies 
( 5.35|) . Apparently we are not sufficiently close to equilibrium, in particular in the low- momentum 
region. Figure 10 shows the evolution of the compound occupation numbers at later times. At NLO, 
redistribution is clearly taking place, and power is being moved towards higher momentum modes. 
In the Hartree approximation nothing much happens. The classical case also shows redistribution 
over the modes and mimics the quantum evolution nicely at low momenta for not too long times. At 
the latest time some concaveness appears to show up, perhaps indicating classical equipartition (cf. 
figured). We come back to this in more detail in the next section. In figure |ll| we show a close-up of 
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Figure 7: As figure || comparing the classical approximation (green/grey) to the NLO (black) result. 
Classically, only the unstable modes are initialized and a continuous function is restored at times \it ~ 5. 
After the transition, the two approximations agree very well. 



the NLO case at the latest time. We include a fit to the low momentum region (neglecting the zero 
mode). As discussed in section [T4] (see figure ||), such a straight line represents a thermalized state, 
with the slope equal to 1/T, and intercept equal to the effective (Goldstone) chemical potential. For 
this fit we find T ~ f.f fj, and fj,^ h ~ 0.32 /i. The temperature is higher than our free-field estimate 
fl3.f| ). This is not surprising since the high- momentum range is still under-populated at this time, 
compared to the thermal state. Figure [12], right, shows the evolution of this temperature in time. 
An exponential fit suggests an asymptotic final temperature of around T/fi= 1.04. 



7.3 Choice of classical initial conditions 

As an aside we comment on the choice of initial conditions for classical simulations of tachyonic 
preheating. Two approaches have been put forward, one used by many authors (see e.g. M, nffl), 



in which all modes are initialized with the "quantum half", i.e. also k > /i in eq. (6.1), the other 
|59| , |f6| , |l4| , |l5| , [ljj, in which only the unstable modes are initialized, as we did in the previous 
section. 

It is clear that initializing all the modes, we risk running into an ultraviolet problem in the 
continuum limit a s — + 0. In the classical approximation, equipartition may cause the divergent 
energy sitting in the high momentum modes to flow into the low-momentum range and generate a 
too high temperature, see []59|], and [^0| for a study with massless gauge fields. However, this did 
not happen on short to intermediate time scales in our 5'f/(2)-rIiggs study in which the gauge fields 
were massive |l^] , and also for scalar fields with only self interactions the problem shows up on a 
much larger time scale than for massless gauge fields [|o) . Thermalization for scalar fields is rather 
slow and it is possible that the high-momentum modes never come into play on time-scales relevant 
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Figure 8: The dispersion relation (cu^/fi 2 vs k 2 /fi 2 ) in the NLO (green/grey) case compared to the 
quadratic approximation (black). 



to numerical simulations of phase transitions. The important criterion must be which initialization 
method reproduces the quantum evolution of the low momentum modes more closely, in terms 
of occupation numbers and equilibration times. We shall take the NLO case as representing the 
quantum result. 

It is also important to keep in mind that, when attempting to calculate topological quantities 
such as the Chern-Simons number in models which include gauge fields, the lattice implementation 
is very sensitive to ultraviolet fluctuations. In simulations of baryogenesis, Chern-Simons number is 
a crucial observable, since it is simply related to the generated baryon asymmetry. This was in fact 
one of the motivations for introducing the initialization restricted to the unstable modes |fl6| , |l4| . 

When initializing all modes, the effective mass receives a contribution from all these modes 
(divergent in the continuum limit), compared to when only the unstable modes are initialized, in a 
way analogous to the quantum vacuum case. We should therefore introduce a counterterm for the 
mass in a way similar to the quantum case Q , 

- -M 2 - A o-^pCi(^i), (7.1) 



where again c\ is given by eq. ( 5.28) ) . For consistency, we also define the occupation numbers 



according to the quantum equation, 



n k + 1/2 = iV-y«^ k )(4^ k ), (7.2) 



rather than through eq. (6.8), and the energy and pressure also have to be renormalized. 

The result is shown in figure [l3|, where the two classical approaches are compared to the NLO 
quantum case. Just after the transition (32 3 , Aq = 1), the classical with only the unstable modes 



- 20 - 




5 10 15 20 25 5 10 15 20 25 5 10 15 20 25 




5 10 15 20 25 5 10 15 20 25 5 10 15 20 25 




Figure 9: The evolution of the dispersion relation (ujI/^j 2 vs fc 2 //i 2 ) in time for NLO, classical and Hartree 
(top to bottom). Overlaid a straight line with slope 1 and intercept (free massless propagator). 



initialized follows the quantum case nicely in the low momentum region. This is also the case when 
initializing all the modes. We conclude that as long as the simple renormalization procedure eq. 
(7.1) is imposed, the number of modes initialized is not crucial at these early times, when classical 
equilibration has not set in yet. 

For longer time (figure [li] left plot, 16 3 , Ao = 6, fj,t = 100), we see that both types of initial- 
ization approximate the quantum behavior, with the 'all initialized' doing somewhat better than 
the 'unstable only'. We do not expect the classical modes with occupation numbers less than one 
to be correct. Surprisingly, the very lowest modes with 3> 1 are also off in both cases. At 
even longer times (figure fli] right plot, [it — 700), the agreement is only qualitative, with either 
intitialization scheme performing badly. For only the unstable modes initialized, we see that the 
occupation numbers are lower than for NLO, presumably because of the lack of scattering with the 
high momentum (stable) modes. In the quantum case there is apparently still some scattering going 
on with the quantum fluctuations of the initial free 'vacuum'. When initializing all the modes, this 
effect appears to be incorporated in the 'all initialized', but it seems to be too large. This could 
also be the reason that the zero mode is lower than for NLO. Eventually, we expect no to approach 
the equilibrium value following from f[||, jslsaj ), n w (v 2 L 3 (n% + l/2)m/N 2 ) 1/2 - 1/2 ~ 19, for 
the parameters of the simulation. 

There is an interesting aspect to this. In figure ^ we see the high momentum modes at the 
same time as figure [14], right; shown is vs k/ fi on a linear scale. In the quantum case, the quantity 
(rik + 1/2) — 1/2 is always positive and here close to zero, as we would expect. This is a consequence 
of the standard quantum commutation relation. The power in the vacuum fluctuations, the "1/2", 
cannot be extracted and moved to other modes. For the classical case with only unstable modes 
initialized, the nu are also manifestly positive numbers, and close to zero. In the classical case with 
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Figure 10: The evolution of ln(l + 1/rik) vs. L0k/[i in time. From top to bottom NLO, classical and 
Hartree. 




kV <VH 



Figure 11: At the latest time pi = 700, the NLO compound dispersion relation and particle numbers 
begin to show the expected behaviour, as in figure ^[ The overlaid lines are (left) asymptotes as in figure]^ 
and (right) a straight-line fit as described in the main text. 



all modes initialized, there is no principle which enforces that (nk + 1/2) should be larger than 1/2, 
and we see that it ends up negative for fc > 3.5 (i. This illustrates the problems associated with 
renormalization in classical field theory out of equilibrium. Only for short times, as we have seen, 
or in equilibrium fell «32|, 63 , such a procedure has a limited validity, albeit with initial-state or 



temperature-dependent counterterms. 
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Figure 12: Left:The evolution of the quantity v //j, as denned in ( 5.37 ). Overlaid the zero temperature 
value (|2.2j). Right: The effective temperature in the low momentum range (from linear fits as in figure [Tl] , 
right). Overlaid, an exponential fit. 




Figure 13: Comparing the NLO result to the classical with unstable modes and all modes initialized. 
fit — 14, just after the tachyonic transition (L/x = 22.4, Ao = 1). 



7.4 Effective temperature in the unstable modes 

In J^, [k| the baryon asymmetry in scenarios based on electroweak preheating was estimated 
using equilibrium concepts, such as the sphaleron rate at an effective temperature and an effective 
chemical potential for Chern-Simons number. The idea was that the large occupation numbers in the 
unstable modes would lead to early equilibration amongst these modes at an effective temperature 
T c g. Such an approach might enable one to bypass the more difficult numerical simulations using 
CP violation in the equations of motion (we note that the comparison made in [ pj did not seem 
to support this idea at a quantitative level). It is therefore of interest to estimate such an effective 
temperature here. 

Equilibrated modes with high occupation numbers have an effective temperature given by the 
classical formula T c g — nk^k, independent of k, and one estimate of the effective temperature 
follows from energy conservation and dominance of the unstable modes |nj , 

r d 3 k jyr cff 37vy 97r 2 9n 2 

7V A <M (2^ w = -6^ = y ° = ^A- ^ Tcff = — » = vm v - (7 - 3) 
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Figure 14: Comparing the NLO result to the classical with unstable modes and all modes initialized. Left: 
fit — 100, after some equilibration. Right: fit = 700. (Lfi = 11.2, Ao = 6). 
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Figure 15: The occupation numbers in the high momentum modes for NLO, classical with unstable 
initialized and classical, all initialized, fit — 700. For the classical, all, n k < for the highest modes, which 
means that n k + 1/2 < 1/2. 



With v = 246 GeV, N = 4, this gives a high T eff ~ 4460 GeV (1820 GeV), for A = 1 (6). The 
actual nkWk of the individual modes as obtained from the (dt4>dt4>) correlator (1/iV Higgs- plus 
(N — 1)/N Goldstone-contributions, as in (5.33)), is shown in figure [l6| (left) at time tfi = 14, for 
the 32 3 lattice with A = 1. We see that the unstable modes have not equilibrated much yet, since 
the effective temperatures vary significantly from mode to mode (the zero mode should be ignored 
as it contains the condensate). Note that the corresponding particle distribution is shown in the 
lower-right panel of figure 0. As the panel tfi = 2.8 in the latter figure shows, the modes near the 
border of the unstable region (k « fi) are not initially highly populated, and also at time tfi = 14 
their effective temperature is much lower than the simple estimate (7.3). The effective temperature 
in the modes near k/fi — 0.5 in figure |l6| (left) is about 100 /i; the average over all unstable (non- 
zero) modes gives T cS = £ k fc<p n k ui k / £ k ,fc <A1 



54/x, about half the estimate (7.3) for N = 4, 
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A = 1. The effective temperatures at later times are shown in figure [16] (right). Equilibration is 
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Figure 16: Left: effective temperatures T c fr,fc = ntuJk (in units of n) of the unstable modes of the 32 2 
lattice at time t[i = 14; the coupling A = 1. Right: same for the 16 3 lattice (still with A = 1), as a function 
of time. Top to bottom: zero mode, L|k|/(27r) = 1, y/2, 

evidently slow even in the unstable momentum region, which is perhaps to be expected because 
of the suppressed scattering of Goldstone bosons at low momenta. In a more realistic setting, the 
gauge fields and other degrees of freedom of the SM will lead to more efficient thermalization and 
lower effective temperatures. In the study [|l2| the SU(2) gauge fields were taken into account, 
which led to effective temperatures T e g « 0.25 — 0.4 ttih « 0.4 — 0.6 \i at trciH ~ 30 — 40 (cf. figure. 
15 and Table 2 in rcf. (l2)), for mu — \/2mw- 



8. Conclusion 

The ^-derived 1 /iV-NLO approximation performed well in this study of a tachyonic phase transi- 
tion, for the 0(4) model at reasonable couplings 9 . It is reassuring that the classical approximation 
agrees so well with the quantum case, not only during the transition, were we used a smaller cou- 
pling (with the accompanying larger particle numbers), but also subsequently for not too long times, 
where we used the larger coupling. In a sense, the classical and NLO approximations support each 
other, the first taking into account all non-linearities without further truncation and the second 
including quantum scattering right from the beginning. For longer times the quantum computation 
is of course superior as it does not suffer from Rayleigh- Jeans ambiguities. 10 

Because of the O(N) symmetry of the initial state, the mean field was zero at all times. Con- 
sequently there was no conflict with Goldstone's theorem - we consider the numerical evidence for 
massless particles to be quite convincing. 

Critical slowing-down of thermalization is expected to take place as a result of the suppressed 
interactions of Goldstone bosons at low momenta. We do see at late times that the intermediate 
range of momenta approach a thermal form, while the very lowest ones may be suffering from 
this critical slowing down. A larger volume (providing very low momentum modes) is needed to 
firmly establish this. In finite volume one would anyway expect a non-zero mean field to vanish as 
the system approaches equilibrium. In this explorative study we did not pursue the finite-volume 
aspects in great detail and more can be done in a finite-size scaling analysis. 

9 At stronger couplings A — > 10 we did see annoying instabilities similar to what we found in approximations based 
on the weak-coupling expansion of <E>. In this sense the larger coupling A = 6 used in Sect. M is not very weak. 
10 For some problems one may have to face a non-uniformity of large N versus large t [p4L 
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